BESSEL SEQUENCES OF EXPONENTIALS ON FRACTAL MEASURES 
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Abstract. Jorgensen and Pedersen have proven that a certain fractal measure v has no infinite set 
of complex exponentials which form an orthonormal set in L 2 (y). We prove that any fractal measure 
fx obtained from an affine iterated function system possesses a sequence of complex exponentials 
which forms a Riesz basic sequence, or more generally a Bessel sequence, in L 2 ()j,) such that the 
frequencies have positive Beurling dimension. 
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1. Introduction 

In |JP98bj . Jorgensen and Pedersen prove two surprising results: 

(i) there exists a singular Borel probability measure a such that there exists a sequence 
{A n }^L C R such that the functions e\ n (x) := e 2mXnX is an orthonormal basis for L 2 (a); 

(ii) The Hausdorff measure v on the middle third Cantor set has the following property: for 
any three {Ai, A2, A3} C R, the set {e\ 1 , e A2 , e A3 } C L 2 {u) is not orthogonal. 

In both cases, the measure arises as the (unique) invariant measure under an iterated function 
system |Hut81j . We prove, in contradistinction to item (ii) above, that every measure \i arising 
from a suitable iterated function system possesses an infinite sequence {X^^Lo such that the 
sequence {eA n }^L * s a Ri esz basic sequence in L 2 (fi). Moreover, this sequence has positive Beurling 
dimension. 

Frames were introduced by Duffin and Schaeffer [DS52] in the context of nonharmonic Fourier 
series, and today they have applications in a wide range of areas. Frames provide robust, basis-like 
representations of vectors in a Hilbert space. The potential redundancy of frames often allows them 
to be more easily constructible than bases, and to possess better properties than those that are 
achievable using bases. For example, redundant frames offer more resilience to the effects of noise or 
to erasures of frame elements than bases. Following Duffin and Schaeffer a Fourier frame or frame 
of exponentials is a frame of the form {e 2mXx }\ & \ for the Hilbert space L 2 [0, 1]. Fourier frames are 
also closely connected with sampling sequences or complete interpolating sequences [OCS02] . 
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Definition 1.1. A sequence {xn} 1 ^^ in a Hilbert space (with inner product (•, •)) is Bessel if there 
exists a positive constant B such that 



X>,x n )| 2 < B\\vf. 



n=l 



This is equivalent to the existence of a positive constant D such that for every finite sequence 
{ci, . . . , ck} of complex numbers 



A" 

E 

n=l 



CnXr 



< D 



A' 



n=l 



Here D 2 = B is called the Bessel bound. 

The sequence is a frame if in addition to being a Bessel sequence there exists a positive constant 
A such that 



#f<El(«^)l 2 <5[HI 2 - 



n=l 



In this case, A and B are called the lower and upper frame bounds, respectively. 

The sequence is a Riesz basic sequence if in addition to being a Bessel sequence there exists a 
positive constant C such that for every finite sequence {ci, . . . , ck} of complex numbers 



C\ 



A 



\ 71=1 



A 

E 

71=1 



< D. 



A 



\ 71=1 



Here C and D are called the lower and upper basis bounds, respectively. 

The main result of Duffin and Schaeffer is a sufficient density condition for {e 27n, ^' x }\^\ to be a 
frame for L 2 [0, 1]. Landau [Lan67j, Jaffard [Jaf91] and Seip [Sei95j "almost" characterize the frame 
properties of {e 2mX x }AeA i n terms of lower Beurling density: 

#(AD[x - h,x + h]) 



V-(A) 



liminf inf 



2h 



Theorem 1.2. For {e 2m> "' x }\ & \ to be a frame for L 2 [0, 1], it is necessary that A is relatively 
separated and P~(A) > 1, and it is sufficient that A is relatively separated and T>~(A) > 1. 

The property of relative separation is equivalent to the condition that the upper Beurling density 

#(An [x-h,x + h]) 



T> + (A) := limsupsup- 

h—^oo xGR 



2h 



is finite. 

For the critical case when T>~(A) = 1, the complete characterization was beautifully formulated 
by Joaquim Ortega-Cerda and Kristian Seip in [OCS02] where the key step was to connect the 
problem to de Branges' theory of Hilbert spaces of entire functions, and this new characterization 
lead to applications in a classical inequality of H. Landau and an approximation problem for 
subharmonic functions. 

In recent years there has been a wide range of interests in expanding the classical Fourier analysis 
to fractal or more general probability measures [DHS091 IHL081 UPOOl IJP98cl IJKS071 ILW061 IL1071 
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IStrOOl IStr06l IYua08j . One of the central themes of this area of research involves constructive and 
computational bases in L 2 (/i), where fj, is a measure which is determined by some self-similarity 
property. These include classical Fourier bases, as well as wavelet and frame constructions. 

For L 2 [0, 1], a sequence of exponentials is Bessel if the frequency set A has finite upper Beurling 
density. For a singular measure v , a necessary (but not sufficient) condition for such a sequence to 
be Bessel in L 2 {v) is that the upper Beurling density of A is (see |DHSW10] ). Since the measures 
we consider here are singular, we shall use Beurling dimension as a replacement for Beurling density. 

Definition 1.3. [CKS08j Let A be a discrete subset of R d . For r > 0, the upper Beurling density 
corresponding to r (or r -Beurling density) is defined by 

V^{A) := limsup sup — . 

h-+oo xeR d ^ 

The upper Beurling dimension (or simply the Beurling dimension) is defined by 

dim + (A) := sup{r > : 23+ (A) > 0} = inf{r > : 23+(A) < oo}. 

Given a set of exponential functions E(A) := {e^ : A G A} we also say that 23+ (A) is the 
r-Beurling density of E(A). 

Definition 1.4. Let R be a d x d expansive integer matrix, B C Z rf , with #B = N > 2. Define 
the iterated function system 

T b (x) = R- 1 (x + b), (x € M. d ). 
For convenience, we let S := R T . 

Let (j>b)beB be a finite set of probabilities, i.e., < p b < 1, J2beBPb = 1- Define the following 
operator T on Borel probability measures on M. d 

(1.1) {T 1 )(E) = Y,P b l(r b -\E)l 

beB 

for all Borel sets E. Equivalently the measure Tj is defined by 

(1.2) / fdT 1 = Y,P^ [ fondj, 

J beB J 

for all continuous functions / on if 

We denote by \i := fiB,p the unique invariant measure for the operator T, i.e. TfJ>B,p = ^B,p-, 
whose existence is guaranteed by [Hu t81j . 

Definition 1.5. For a Borel probability measure 7, if A = {A n }~ =0 C M is such that {e A J C L 2 ( 7 ) 
is a Bessel sequence, we say A is a Bessel spectrum for 7. Likewise, A is a Riesz basic spectrum if 
{e\ n } is a Riesz basic sequence in L 2 {^f). 

In the classical Lebesgue measure case, it is relatively easy (with the help of Theorem 1 1 . 2 [) 
to construct frames/Riesz bases or more generally Bessel sequences/Riesz sequences {e 2mX ' x }AeA 
with A having positive Burling density. However, this is not the case anymore for fractal measures. 
Indeed, for the fractal measure Hb& in the case that R = 3, B = {0,2}, and Po = P2 = 1/2, the 
corresponding measure ^3 (which is the Hausdorff measure on the middle third Cantor) has the 
property that {3™ : n = 0, 1, . . . } is NOT a Bessel spectrum (and hence can not be a Riesz basic 
spectrum) [DHSW101 Proposition 3.10]. Note that this set {3 n : n = 0, 1, . . . } is very "sparse" and 
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in fact it has the Beurling dimension equal to 0. One of the open problems for the fractal measure 
/i3 is that whether frames or Riesz bases spectrum exist. In [DHSW10, Theorem 3.5] it was proved 
that for a fractal measure fJ>B pi a necessary condition for A to be a Bessel spectrum is that the 
Beurling dimension of A is at most log^ B, and that the Beurling dimension of A is equal to log R B 
(under a mild technical condition) in order for A to be a frame spectrum. The above example 
(A = {3 n : n = 0,1, . . .} ) shows that this finite Beurling dimension condition is not sufficient for 
A to be even a Bessel spectrum. This naturally leads to the existence problem for Bessel spectrum 
and Riesz basic spectrum with positive Beurling dimensions. The main purpose of this paper is to 
prove that Bessel spectrum and Riesz basic spectrum with positive Beurling dimension exists for 
all the fractal measures lib, p - We believe that this is an important positive step toward answering 
the question of whether frames or Riesz bases spectra exist for the fractal measure 113. 

2. Spectra of Positive Beurling Dimension 

We start with our main theorem: 

Theorem 2.1. Let R be adxd expansive integer matrix, S B C 7L d , (pb)beB a list of probabilities 
and let ll = lib, p be the invariant measure associated to the iterated function system 

Tb (x) = R- l {x + b), (xeR d ,beB) 

and the probabilities (pb)beB- Then ll has an infinite Riesz basic spectrum of positive Beurling 
dimension. 

The proof proceeds via a series of lemmas. Throughout the remainder of the paper, R, B, p are 
fixed, and /i := lib, p - 

Lemma 2.2. The Fourier transform of the invariant measure fi satisfies the scaling equation 

(2.1) Jl(x) = miS^x^iS^x), (x E R d ) 
where 

(2.2) m(x) := ^p b e 2 ™ b x , (x € R d ) 

beb 

The function fx is given by the infinite product formula 

00 

(2.3) fi{x) = Y[m ^S- k x^j , (x € M. d ) 

fc=i 

The infinite product converges uniformly on compact subsets. 

Proof. Apply the Fourier transform to the invariance equation (II. 2D . See e.g. |JP98ai IDJ06] for 
details. □ 

Lemma 2.3. There exists p G N, < p < 1, and a finite set A C Z d \ {0} with ffA > 2 such that 
ifM:=max{\\S- p a-S- p a'\\:a,a'eAu{0},a^a'},then 

(2.4) \m(S~ p {a-a') + x)\< p 

for all x with \\x\\ < ^pprj^> an d f or oil a,fl' € AU {0} with a 7^ a 1 . 
In addition, the elements of A are incongruent modS p Z d . 
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Proof. We have to stay away from the points where \m\ is 1. Note that |m(aj) | = 1 implies 



beB 



1. 



Since G B, the term po ■ 1 appears in the sum. Using the triangle inequality we must have 
e 2mb-x _ j f or a n 6 g -B and therefore 6 • x G Z. It follows that 

5i := {x : |m(x)| = 1} = {x : b ■ x G Z}. 

Since «Si has Lebesgue measure zero in M. d , we can find two distinct points xq,x\ such that 

ll^oll) \\ x i II ^ 1 with ±^0; i^ij i(^o — %i) & S\. 

Choose 5' < 1 and let c := 1 1 jS* 1 1 1 < 1 (since the matrix S is expansive). We can pick p G N large 
enough such that 

A r P 

(2.5) -p— < 5'. 

y ' 1 - CP 

Moreover, since the volume of the lattice S~ p 7L d goes to as p gets large, we may choose p so that 
additionally there exist integers Oq ^ ai G Z d with ||z„ - 5" " p ao||, ||xi — S p a\\\ < 5'. 

Let A := {oo,ai}. Then ||S _p ai|| < 1 + 5' < 2 so M as defined in the hypothesis will be less than 
4. If \\y\\ < Mc p /(1 — c p ) then for a, a' G vl U {0}, a 7^ a', there exists x,x' G {xo,xi,0} such that 

\\(S- p (a-a')+y) - (x - x')\\ < \\S' p a - x\\ + \\S~ p a' - x'\\ + \\y\\ <35'. 

Thus, if 5' is small enough, S~ p (a — a') + y, being close to x — x', stays away from the set Si so by 
uniform continuity of m, there is a p < 1 such that 

\ m (S- p (a-a') + y)\ < p, 

for all y with \\y\\ < Mc p /(1 — c p ). This proves the existence of p, p and A. 

The elements in A U {0} cannot be congruent modS' p Z rf because |m(S' _p (a — a'))\ < 1, while 
m(k) = 1 for k G Z rf and m is Z rf periodic. 

□ 

Definition 2.4. Let p, p and A be as in Lemma 12 .31 Let 

A(A,p) := {a + S p ai + ■■■+ S pr a r : 0i G A U {0}}. 

We identify the integer ao + 5" p ai + • • • + S pr a r with the word a^ai . . . a r and with the infinite word 
a^ai . . . a r a r+ \ . . . , where = for i > r + 1. Since the elements in A U {0} are incongruent 
m.odS p 7 l d , different digits means different integers. 

For two such A = a^ai . . . , A' = a' a'i . . . , we define the Hamming distance between them by 

d p (X,X') = #{i:a l ^a' i }. 

Lemma 2.5. Let p, p, A, M be as in Lemma \2.3[ Let A = a^ai . . . , A' = al^a'^ ... be distinct words 
with digits in A. Then 

\p(X-X')\ </p( A ' A '). 

Proof. We use the infinite product formula for p, and we group every p terms. We have 

00 

p(x) = H m^(S- np x), (x G R d ) 

n=l 
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where 

rrS p \x) = m{x)m{Sx) . . . m(S p ~ 1 x). 
Since \m\ < 1, we get for any I C N, 

(2.6) |£(z)| < l\ \m^{S- np x)\. 

n£l 

Suppose a n / a' n . For k > n, S~ np (S kp (ak — a' k )) G Z d . Therefore, we have 

S~ np {\ - A') = 5- p (a n _i - + S- 2p {a n - 2 - a' n _ 2 ) + ■■■ + S~ np {a - a' ) mod Z d , 

and since 

Mil 9 — 1 11^ 

\\S- p (a n _ 1 -a' n ^)+---+S- np (ao-a' Q )\\ < \\S^\\(M+\\S^\\'M+- ■ ■+\\S^ n - 1 > ||-M) < j _' |S _J [p 
with Lemma 12.31 we obtain 

|m(,S~ np (A- A')) | < p. 

Thus, using (12. 6h . we obtain 

|M(A - A') | < JJ \m^(S- np {\ - A') | < /f( A ' A '), 

where / := {n : a n ^ a^} with = d p (A, A'). □ 

Remark 2.6. By changing the set of digits A and p we can assume that p is as small as we want. 
This is because we can replace each digit in A by a repetition of it, say / times. So, for example 12 
is replaced by 111222, where / = 3. By doing this the distance between any two words is multiplied 
by I. So if we replace p by p ■ I and A by A® := {a"' := aa . . . a : a G A}, the number p in Lemma 

i times 

12.51 is replaced by p l , which can be made as small as needed. More precisely, we have 

<WA« 7 (0 ) = I ■ d pl>Am (\,i), for A, 7 G A« 
The distance , A (i) counts each digit as just one. 

Lemma 2.7. Letp,p,A be as in Lemma \2.'J[ and let A C A{A,p). Suppose 

C:=sup P dp{X ' X>) < oo. 

AeA A'eA\{A} 

Then A is a Bessel spectrum with Bessel bound 1 + C . 
Moreover, if C < 1, then A is a Riesz basic spectrum. 

Proof. Using Lemma 12.5^ we have for all A £ A: 

Im( a - A ')i < 1 + pdp{x ' x,) ^ i+c - 

A'eA A^A' 

An application of Schur's lemma shows that the Grammian of the set {e^ : A G A} is bounded, 
and, if C < 1, it is also diagonally dominant, hence invertible (see |Chr03l Proposition 3.5.4]). 
This implies that A is a Bessel spectrum with bound 1 + C and is a Riesz basic spectrum when 
C<1. □ 
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Lemma 2.8. Let A be an alphabet with 2 letters. Then there exists ko > 1 such that for every 
k > Uq and every n there is set A n containing 2" words of length kn, i.e., A n C A kn , such that the 
Hamming distance between any two distinct words in A n is at least n. 

Proof. This is a consequence of the Gilbert- Varshamov bound [vL99] which states that if A q (m,d) 
is the maximum possible size of a q-aiy code C with length m and minimum Hamming distance d 
(a g-ary code is a code over the field ¥ q with (/-elements), then 

A q (m, d) > 



E?=doU<7-iy' 



where C m are the binomial coefficients. 

In our case, q = 2, d = n and m = kn. Let H(x) 
argument in [BIPW10J we can show that for k > 2 



-x log 2 x 



[1 — x) logofl — x). Using the 



Y,C m <2 H ^ kn . 

j=0 



In fact, let u = j. Then 



-H(u) 



and so 



i = (u+(i-u))*»>x;c*y(i 

3=0 



1 - U 

- u) kn ~ j 



>E C L 

3=0 

Thus we get 



1-u 



:i- u 



,kn 



3=0 



l-U 



(1-u) 



i=o 
(1-u) 



ii 



1 - n 

kn 



(1-u) 



kn 



3=0 



3=0 



hence 



Since 1 - H(l) 



->kn 



A2(kn, n) > 



> 



7=0 u fa Z^i=0 ^fc 



> 2 ( 1 -^(i)) fc «, 



1 as A; — > oo, there exists a feo such that for k > ko, (1 — H(^))k > 1 and 
A 2 (kn,n) > 2 n . 



□ 



Proof of Theorem \2.1\ To complete the proof of Theorem 12.11 we construct a set A that satisfies 
the hypothesis of Lemma 12.71 

Let A, p, p as in Lemma 12.31 where A has two non-zero elements. Without loss of generality, by 
Remark 12.61 we may assume that p < 1/4. 

We choose a sequence qi,q%,... of natural numbers such that q\ -\ + q n -\ + 1 < q n , for all n; 

for example q n = 2 n . 

Let k > ko where ko is as in Lemma \2.8\ let A n be the set of words of A kqn guaranteed by Lemma 
with at least 2 qn elements and the Hamming distance between any two words is at least q n . 
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We define A by concatenating words in Aj as follows 

A := {Ai . . . A n : Aj e k h n > 1}. 
Fix A = Ai . . . Xn; we wish to estimate 

a'ga\{a} 

For any A' = A^ . . . X' m £ A, there exists a natural number r(A') which is the largest index such 
that A r ^ X' r . The Hamming distance between A r and A', is at least <2V(A')> so the Hamming distance 
between A and A' is also at least q r (\')- Thus, for a fixed ro € N, we count how many A' € A which 
have r(X') = vq. If ro < N, the number of possibilities is 

giH — h?r 

2? 1 2 91+,?2 -)-... + 2 9lH — '" 9r o < ^2 2* < 2 91 " 1 — ^ro+i < 2 2l?r o = 4 9r o . 

i=0 

If ro > N, then A n+ i = • • • = A ro = = A' +1 with A^. ^ 0, so the number of possibilities is at 
most 

It follows that 



2<ll H \~<Ir < 2 2lJr o = 4? r o 



p<h{\x) < A qr p qr . 

A'GA\{A} r=l 

Since p < 1/4, this sum converges. 

We now associate words in A with integers written in base S p with coefficients from the words 
in A as in Definition 12.41 Combining Lemmas 12.51 and 12.71 we conclude that A is a Bessel spectrum. 
Taking qi larger if needed, we can get the sum to be less than 1, so we obtain a Riesz basic spectrum. 

It remains to prove that the Beurling dimension is positive. We will use the following lemma 
that can be obtained by a straightforward computation. 

Lemma 2.9. To compute the Beurling dimension, the unit cube [—l,l} d in Definition \1.3\ can be 
replaced by any bounded set Q that contains in the interior. 

Consider the elements in A\ . . . A n . There are at least 2 Qn such elements, since just A n has 2 9n 
elements. The length of such a word is kqi + • • • + kq n < 2kq n — 1. Let C = max a ^A Then the 
integer represented by this word will have absolute value less than 

C + \\S\\ P C + ■■■ + \\S\\ P ' (2*9"-i)c < DUS'll 2 ^ 

for some constant D. 

Therefore, in the ball of radius D||S'|| 2A: ' ?nP there are at least 2 q " elements in A. Using Q = B(0, D) 
in Lemma 12.91 x = 0, and h = \\S\\ 2kqnP in the Definition 11.31 this implies that the Beurling 
dimension of A is at least logngpfcp 2 > 0. □ 

We conclude with a remark concerning the usual Cantor middle third set. This set, and its 
invariant measure, is generated by the iterated function system with parameters R = 3, B = {0, 2}, 
and po = p2 = The invariant measure ^b,p f° r these parameters is the measure v mentioned 
in item {ii) in the introduction; for this measure there are no three pairwise orthogonal complex 
exponentials, and hence this measure possesses no orthonormal sequence of exponentials [JP98b]. 
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However, by applying Theorem 12.11 this measure on the Cantor set does possess a Riesz basic 
sequence of complex exponentials. 
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